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Recently, magnons, which are quasiparticles describing the collective motion of spins, were found to 
undergo Bose-Einstein condensation (BEC) at room temperature in films of Yttrium Iron Garnet (YIG). 
Unlike other quasiparticle BEC systems, this system has a spectrum with two degenerate minima, which 
makes it possible for the system to have two condensates in momentum space. Recent Brillouin Light 
Scattering studies for a microwave-pumped YIG film of thickness d= 5 jum and field H = 1 kOe find a 
low-contrast interference pattern at the characteristic wavevector Q of the magnon energy minimum. In this 
report, we show that this modulation pattern can be quantitatively explained as due to unequal but coherent 
Bose-Einstein condensation of magnons into the two energy minima. Our theory predicts a transition from 
a high-contrast symmetric state to a low-contrast non-symmetric state on varying the d and H, and a new 
type of collective oscillation. 

Bose-Einstein condensation (BEC), one of the most intriguing macroscopic quantum phenomena, has been 
observed in equilibrium systems of Bose atoms, like 4 He 1,2 , 87 Rb 3 and 23 Na 4 . Recent experiments have 
extended the concept of BEC to non- equilibrium systems consisting of photons 5 and of quasiparticles, 
such as excitons 6 , polaritons 7 " 9 and magnons 1011 . Among these, BEC of magnons in films of Yttrium Iron Garnet 
(YIG), discovered by the group of Demokritov 11 " 17 , is distinguished from other quasiparticle BEC systems by its 
room temperature transition and two-dimensional anisotropic properties. In particular, the spin-wave energy 
spectrum of a YIG film shows two energetically degenerate minima. Therefore it is possible that the system may 
have two condensates in momentum space 18 . An experiment by Nowik-Boltyk et al 17 indeed shows a low- 
contrast spatial modulation pattern, indicating that there is interference between the two condensates. Current 
theories 19-24 do not describe the appearance of coherence or the distribution of the two condensates. 

This report points out that a complete description of BEC in microwave-pumped YIG films must account for 
the 4th order interactions, including previously neglected magnon-non-conserving terms originating in the 
dipolar interactions. The theory explains not only the appearance of coherence but also quantitatively explains 
the low contrast of the experimentally observed interference pattern. Moreover, the theory predicts that, on 
increasing the film thickness from a small value d, there is a transition from a high- contrast symmetric state for 
d < d c , with equal numbers of condensed magnons filling the two minimum states, to low-contrast coherent non- 
symmetric state for d > d c , with different numbers of condensed magnons filling the two minimum states. In 
comparatively thin films (d < 0.2 fim) the same transition can be driven by an external magnetic field H. At 
another critical thickness d* > d a the sum of phases of the two condensates changes from n to 0; at this transition 
point the system is in a completely non -symmetric state with only one condensate, for which there is no 
interference. In the experiment of Ref. 17 the thickness of film was larger than d*. We suggest that the phase 
transitions may be identified by measuring the contrast of the spatial interference pattern for various d and H. We 
also predict a new type of collective magnetic oscillation in this system and discuss the possibility of domain walls 
in non- symmetric states. 

Results 

Phase diagram. We consider a YIG film of thickness d with in-plane magnetic field H (see inset of Fig. 1 ) . The 4-th 
order interaction of condensate amplitudes reads 25 " 27 : 

V 4 = A^c^cqCq + c^qC^qC-qC-q^ +2Bc^_ Q C- Q c Q + Cy [ Q CQCQC-Q + J_ Q C-QC-QCQ + hx^. (1) 

where, c ±Q and c\ Q are the annihilation and creation operators for magnons in the two condensates located at the 
two energy minima (0, ±Q) in the 2-D momentum space (see. Fig. 1). The coefficients in Eq.(l) are: 
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Figure 1 | The magnon spectrum in the k z direction for d = 5 \im and 
H = 1 kOe. The inset is a schematic diagram of YIG film. 
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with oc 1 = u 4 + 4u 2 v 2 + v 4 , a 2 = 2wV and a 3 = 3uv(u 2 + v 2 ). Here, w 
and v are the coefficients of Bogoliubov transformation (see the 
Methods section for details). S = 14.3 is the effective spin, N the 
total number of spins in the film, M the magnetization and 
HwM = y4nM with gyromagnetic ratio y = 1.2 X 10~ 5 eV/kOe. D 
is proportional to the exchange constant and F k = (1 — e~ kd )lkd. 
Similar results for the coefficients A and B were obtained in Ref. 19. 
Coefficient C, which violates magnon number conservation, was not 
considered previously. Below we show that C is the only source of 
coherence between the two condensates. The three coefficients A, B 
and C, whose values as functions of H are shown in Fig. 2 for two 
typical values of d, determine the distribution of condensed magnons 
in the two degenerate minima. Ref. 19 assumed a symmetric state 
with condensed magnons in both minima having equal amplitudes 
and equal phases. Later, Ref. 20 assumed filling of only one 
minimum. More recently Ref. 24 considered Josephson-like 
oscillations by starting from two condensates with equal numbers 
of magnons but different phases. Our theory predicts coherent 
condensates and the ratio of their amplitudes without any 
additional assumptions. 

In terms of condensate numbers N ±Q and phases <j> ±i the con- 
densate amplitudes are c+q = ^/N+Qe 1 ^ . Substituting them into 
eq.(l) we find: 



=a(n 2 q +nI 

+ 2CcosO>( / 



+ 2BN Q N_ Q 



N J Q Nl 



(3) 



where we introduce the total phase O = </>+ + </>_. To minimize this 
energy, O must equal n for C> 0 and must equal 0 for C < 0. Fig. 2 
shows that the sign of C changes for different d and H, which indi- 
cates a transition of <I> between 0 and n. For both C> 0 and C < 0 a 
coherence between the two condensate amplitudes is established. In 
contrast to the Josephson-like interaction, the sum rather than the 
difference of the two condensate phases is fixed. 
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Figure 2 | The interaction coefficients A, B and C (in units of mK/N, 
with Nthe total number of spins in the film) as a function of magnetic 
field H for film thickness (a) d = 1.0 um and (b) d = 0.1 um. 

Since the total number of condensed magnons N c = N Q +N- Q is 
uniquely determined by the pumping (see Methods), the energy is 
minimized only by the so far unspecified difference 3 = N Q — N- Q . 
In terms of N c and 3 the condensate energy eq.(3) is: 



Va 



(A + B)N 2 -(B-A)3 z -2\C\N c JN 2 -3 



(4) 



The ground state of the condensates depends on the criterion 
parameter: 

A = A-B+\C\. (5) 

When A > 0, 3 = 0 minimizes the energy, so the two minima are 
filled with equal number of condensed magnons. This is the sym- 
metric state with N Q = N- Q . When A < 0, the minimum shifts to 

3 2 C 2 

— - = 1 5-. This is the non- symmetric state with N Q ^ 

N 2 (B — A) 7 Q 

N-q. The transition from symmetric to non-symmetric state at 
A = 0 is of the second order. There is no metastable state of these 
states. At C = 0 one finds 3 = ±N C , which corresponds to a com- 
pletely non- symmetric state with only one condensate. The ground 
state of the non- symmetric state is doubly- degenerate, correspond- 
ing to the two possible signs for 3. Fig. 3 shows that for a film 
thickness of about 0.05 ^m, the symmetric state is energy favorable 
up to H = 1.2 T. For d = 0.08 ^m, on increasing H to about 0.6 kOe, 
there is a transition from symmetric to non-symmetric state. For a 
larger thickness d = 0.1 ^m or d = 1 ^m, the ground state is non- 
symmetric for H > 0.3 kOe. 

Fig. 4 shows that the phase diagram in (d, H) space has three 
different regions, separated by two critical transition lines, d c (H) 
and d*(H), corresponding to A = 0 and C = 0, respectively. For 
d= 0.1 3— 0.16 ^m, the system possesses re-entrant behavior (NS 0=7i, 
to NS <D=0, to NS 0=7c). As shown below, measurement of the 
contrast, or modulation depth 17 , of the spatial interference pattern 
permits identification of the different condensate states. 

Zero sound. In two -condensate states the relative phase d<j> = <j>+ — 
(/)- is a Goldstone mode. Its oscillation, coupled with the oscillation 
of the number density 3n = n Q — n- Q represents a new type of 
collective excitation, which we call zero sound (as in Landau's 
Fermi liquid, this mode is driven by the self- consistent field rather 
than collisions). Solving a properly modified Gross-Pitaevskii 
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Figure 3 | The criterion of transition from non- symmetric to 
symmetric phase, A (in units of mK/N), as a function of magnetic field H 
for different values of thickness d. 

equation (see Methods), we find its spectrum. In the symmetric state 
it is: 



h 2 k A , k 2 
— +N C A-. 
Am 1 m 



(6) 



The magnon effective mass is of the order of the electron mass. The 
density of condensed magnons n c = NJ Vis about 10 18 cm" 3 and A ~ 
10 mK/N. The sound speed for small k in this case is v 05 = ^/N c A/m, 
which is about 100 m/s. Near the transition point A = 0, the velocity 
of this zero sound goes to zero. For the non-symmetric state, the 
spectrum is: 



ft 2 A: 4 , w ,k 2 
— K +N C (B-A)( K -1)-, 



where k 



_ (B-A) : 
C 2 



. In the experiment of Ref. 17, K ~ 10 4 and B — A 

= 8.4 mK/N. An estimate of the sound speed gives 3 X 10 3 m/s. The 
dispersions of zero sound for symmetric and non-symmetric states 
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Figure 4 | The phase diagram for different values of thickness d and 
magnetic field H. 



are shown in Fig. 5. Note that the range of applicability of the linear 
approximation strongly shrinks at small C since the density of one 
of the condensates becomes very small and the phase fluctuations 



Domain wall. Since the ground state of the non -symmetric state is 
doubly degenerate, it can consist of domains with different signs of 3 
separated by domain walls. The width w of a domain wall is of the 



order of 



h 2 



2mN c \A 
wall energy per unit area is 



. For the data of Ref. 17, w ~ 10 jam. The domain 
10" 9 J/m 2 . 



Discussion 

The ground state wave function *F(z) generally is a superposition of 
two condensate amplitudes ^(z) 
c ±Q = ^N ±Q e i( ' ' 



(c Q e iQz + C- Q e~ iQz ) /y/V, where 
and Vis the volume of the film. The spatial struc- 
ture of ^(z) can be measured by Brillouin Light Scattering (BLS). 
The BLS intensity is proportional to the condensate density 
| x F| 2 = n Q + n_ Q + 2 v /n Q n_ Q cos(2Q2: + ^ + -0_). 

In their recent experiment, Nowik-Boltyk el al 17 observed the 
interference pattern associated with the ground state. They found 
that the contrast of this periodic spatial modulation is far below 
100%; of the order 3%. The present theory can quantitatively explain 
this result. In the experiment of Ref. 17, d = 5.1 fim and H = 1 kOe, 
eq.(2) for A, B and C then gives A = -0.168 mK/N,B = 8.218 mK/N 
and C = —0.203 mK/N, so A < 0. This corresponds to the non- 
symmetric state, where the ratio of the numbers of magnons in the 
y N_o C 2 6 



two condensates is - 

o I \max I I mm 

I^FI 2 +| X F| 2 -' 

m/JY 1 \rnin 



Nn 



4(B-A) : 



(assume 3 > 0). The contrast is 



Since C<3 and iV_Q<iV Q , ^2 



Nn 



r . For the above values of A, B and C, B is of order 2.4%, in 

\B-A\ 

good agreement with experiment. The smallness of C (and A) in 
comparison to B is associated with a large parameter d/l where 

/ = \ / — is an intrinsic length scale of the system and / ~ 

V^ M | C | 

10 6 m. In terms of this parameter, — ' 

Experimentally 17 the contrast p reaches the saturation value at a 
comparatively small pumping power, corresponding to the appear- 
ance of BEC. This agrees with our expression for which depends 
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Figure 5 | Dispersion of zero sound as a function of wave vector in the 
direction of external magnetic field for symmetric and non- symmetric 
states, respectively. For the non-symmetric state, we choose H = 1 kOe 
and d = 5 jam. 
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and NS denote symmetric and non-symmetric state, respectively. 

only on film thickness d and magnetic field H. By varying d and H, 
the contrast can be changed. Specifically, in the symmetric state, 
P = 1; in the non-symmetric state, P < 1; and in the completely 
non-symmetric state with only one condensate, P = 0. Therefore, 
measurement of the contrast for different values of d and H can give 
complete information on the phase diagram of the system, for com- 
parison with the present theory. 

Fig. 6 plots C, A and P as functions of the film thickness d at fixed 
magnetic field H = 1 kOe. For small dasH increases the system is in 
the high-contrast symmetric state. At a larger thickness d c = 
0.07 /im, the sign of A changes, indicating a transition from the 
symmetric to the low-contrast non-symmetric state. As d further 
increases, to d* = 0.17 ^m, C changes sign, where the total phase 
O changes from n to 0. Only at this point d* does the zero-contrast 
state (with only one condensate) appear. Correspondingly, a char- 
acteristic cusp in the contrast P appears near d*. 

To conclude, we have calculated the 4-th order magnon- 
magnon interactions in the condensate of a film of YIG, including 
magnon-non-conserving term responsible for the coherence of two 
condensates. We predict a phase transition from symmetric to 
non-symmetric state that happens at a reasonable magnetic field 
> 0.2 kOe in sufficiently thin YIG films d < 0.1 fim. We also predict 
that within the non-symmetric state there is a thickness d*(H) where 
the modulation in the observed interference pattern should totally 
disappear. 

Methods 

Magnon spectrum. In a YIG film with an in-plane external magnetic field H, the 
magnon dispersion has been studied extensively 28 " 30 . At low energies, YIG can be 
described as a Heisenberg ferromagnet with effective large- spin S = 14. 3 19 24 on a 
cubic lattice. The Hamiltonian consists of three parts: 

n= -jJ2 s i' s i+ H D-ynJ2 sz i> W 

(y) 

the nearest neighbor exchange energy, the dipolar interaction and the Zeeman energy. 
We take y to be perpendicular to the film and the magnetic field to be in the plane 
along z. It is convenient to characterize the exchange interaction by the constant 
D = 2JSa 2 = 0.24 eVA 2 . The dipolar interaction can be calculated using the method 
indicated in Refs. 20,30. The competition between the dipolar interaction and 
exchange interaction leads to a magnon spectrum co k with minima located at the two 
points in 2D wave-vector space given by k = (0, ± Q) (i.e. along z), with an energy gap 
A 0 . For film thickness d = 5 fim and magnetic field H = 1 kOe, we find 
Q = 7.5 X 10 4 cm" 1 and A 0 = 2.7 GHz. In the experiment 17 the wave vector Q was 
found to be about 3.5 X 10 4 cm" 1 , i.e. about half the predicted value. The reason for 



this may be associated with a rather shallow minimum. In such a situation small 
corrections to our approximate formula can have a large effect on the value of Q. The 
lowest band of the magnon spectrum can be calculated with the help of the Holstein- 
Primakoff transformation 32 expressing the spin operator S in terms of boson 
operators a and a f . 

To second order in a and a 1 ", the Hamiltonian eq.(8) is: 



*»=£ 



A k ala k +^B k a k a- k +^Bla\a ] _ i 



with 



A k = yH 0 + Dk 2 + y2nM( l-F k ) sin 2 9 + y2nMF k 
B k = y2nM(l-F k )sm 2 6-y2nMF k 



(9) 



(10) 



where F k = (1 — e~ kd )lkd and M is the magnetization of the material (4nM = 
1.76 kG). Here, 6 is the angle between the 2D wave vector k and the direction of 
magnetic field (z). Ho of eq.(9) is diagonalized by the Bogoliubov transformation 



t - + u (A k + hcD k 
a k = u k c k + v k c ] _ k with u k = I 2 

leading to the magnon spectrum: 



1/2 



and v k = sgn(B k ) 



A k -h(p k 
2hw k 



1/2 



hco k = (A 2 k -\B k \ 



2X1/2 



(11) 



Fig. 1 gives the magnon spectrum along k z for typical values of thickness d and 
magnetic field H. 

Number of condensed magnons. N c = Nq + N-q. Experimentally, the spin lattice 
relaxation time is of order 1 fis, whereas the magnon-magnon thermalization time is 
of order 100 ns; the magnons are long-lived enough to equilibrate before decaying, 
thus making BEC possible 11 . After the thermalization time the pumped magnons go 
to a quasi-equilibrium state with a non-zero chemical potential fi. The number of 

pumped magnons N p = N(T, fi) - N(T, 0), where N(T,fi) = ]T k e{oJk _^ l)/T _ l , is 

determined by the pumping power and the magnon lifetime, ji is a monotonically 
increasing function of N p but cannot exceed the energy gap A 0 . Therefore, on further 
increase of pumping some of the pumped magnons fall into the condensate. The 
equation N pc = N(T, A 0 ) — N(T, 0) thus defines the critical line of condensation. Since 
A 0 <C T and N p -C N(T, 0) this equation can be satisfied at a rather high temperature. 
The total number of condensed particles is 11 31 

N c =N p -N(T,Ao) +N(T,0) =N(2>) -N(T,A 0 ). (12) 

In exactly 2D systems BEC formally does not exist since in the continuum 
approximation the sum in N(T, ft) diverges. However, for strong enough pumping the 
chemical potential approaches exponentially close to the energy gap: A 0 — [i ~ A 0 
exp(-Np/N 0 ), where N 0 = VTm/h 2 . For N p /N 0 > ln(T/A 0 ) all pumped magnons 
occupy only one or two states ± Q. 

Eq.(12) determines only the total number of particles in the condensate. The 
distribution of the condensate particles between the two minima remains 
undetermined in the quadratic approximation. To resolve this issue it was necessary 
to consider the fourth order terms in the Holstein-Primakoff expansion of the 
exchange and dipolar energy. Terms of third order in this expansion occur due to the 
dipolar interaction, but they vanish for the condensate values of momentum (0, ± Q) 
since in the third order the total momentum cannot be zero. 

Zero sound. Here we provide details in calculating the zero sound spectrum. Let us 
consider small deviations from the static symmetric solution n Q = ri- Q = n c /2, 
(f>+ = n — (f>- = 0 so that n ± Q = nJ2 + Sn ± with Sn + = — Sri- = Sn/2 and 
d(j)+ = -§(j)- = 5<j)l2. Then 



W + | 2 + 1 W_ | 2 ) +AV(|*F + | 4 + |*F_ | 4 ) 

+ 25y| x F + | 2 | x F_| 2 + CV( x F + x F_+ x F* + x F*_)(| x F + | 2 + | x F_| 2 
On linearizing, the energy reads: 
E-- 



= \dr(— \S/5n\ 2 + h ^ |V^| 2 + — 5n 
J \4mn c u - 



Using the commutation relation [ 
ihAcf) = [A0,H], we obtain: 

^- 

dt 



4m 2 
Sn] = /, and the equation of motion 

h 2 



2mn c 



-V 2 &n + AVSn, 



(B) 
(14) 

St 2m— (15) 
Taking Fourier transforms of the above two equations in coordinate and time, one 
arrives at dispersion relations Eq.(6). 
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